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Abstract.—The subalpine rain páramo of Isthmian Central America is an area with a high level of endemism. 
The salamanders of the Bolitoglossa subpalmata Species Group are restricted to the highlands of Costa 
Rica and Panama, including the subalpine rain páramo. During explorations of the páramos in the Cordillera 
de Talamanca, we found populations of Bolitoglossa that were referred to B. subpalmata S.G. These new 
populations were compared phylogenetically with all species within the species group using the 16S and cyt b 
mitochondrial genes; and they were also compared morphometrically with B. kamuk and B. pesrubra. Herein 
we described a new species of the B. subpalmata Species Group, which inhabits the subalpine rain páramo and 
montane forest surrounding the páramo in the Southeastern region in Costa Rica. This new species highlights 
the role of the subalpine rain páramo in the speciation of salamanders in the highlands of Isthmian Central 
America, as well as the need to protect this unique and vulnerable habitat. 
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Resumen.—El paramo subalpino de América Central istmica es alto en endemismos. Las salamandras del 
grupo de especies Bolitoglossa subpalmata están restringidas a las tierras altas de Costa Rica y Panamá, 
incluyendo el páramo subalpino. Durante exploraciones a los páramos en la Cordillera de Talamanca 
encontramos varias poblaciones de Bolitoglossa que fueron asignadas al grupo de especies B. subpalmata. 
Estas nuevas poblaciones fueron comparadas filogenéticamente con todas las especies dentro del grupo de 
especies de B. subpalmata; además fueron comparadas morfométricamente contra B. kamuk y B. pesrubra. 
Aquí se describe una nueva especie del grupo de especies B. subpalmata, habitante del páramo subalpino y 
del bosque montano cerca del páramo en la región sureste de Costa Rica. Esta nueva especie resalta el rol 
del paramo subalpino en la especiación de las salamandras de las tierras altas de América Central Ístmica y la 
necesidad de proteger este hábitat ünico y vulnerable. 
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Introduction 


The subalpine rain paramo from Isthmian Central 
America (ICA) is an extremely reduced ecosystem that 
is fragmented and vulnerable. In ICA, the subalpine rain 
paramo is isolated in the highlands, especially in the 
Cordillera de Talamanca’s summits, and covers less than 
150 km? (Kappelle and Horn 2016). The paramo in ICA can 
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be considered a “sky island” located in the northwestern 
part of the Neotropical paramo, separated by straight-line 
distances of 1,100 km to Colombia’s Paramo and 1,000 
km to the highlands of Nuclear Central America. The 
subalpine paramo of ICA is characterized by low species 
diversity but a high level of endemism, especially among 
plants, although several vertebrate species are exclusive 
to the paramo (Kappelle and Horn 2016). Even though 
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Fig. 1. (A) Map showing the known populations of the species within the Bolitoglossa subpalmata species group. (B) Map showing 
the localities for B. bolanosi sp. nov. and the species close to it in the páramo of the Cordillera de Talamanca. The open shapes 
indicate the historical data and the solid shapes indicate the localities with molecular data included herein. 


the subalpine rain paramo of ICA is relatively small in 
area, several regions remain unexplored and several 
species may possibly remain unnamed. To the best of 
our knowledge, only three species of salamanders are 
exclusively found in the ICA’s paramo: Bolitoglossa 
kamuk, B. pesrubra, and B. pygmaea. 

Most salamanders of the ICA are associated with the 
highlands, especially on the Cordillera de Talamanca, 
where studies have documented high species turnover 
between sites and along elevational gradients (García- 
Paris et al. 2000; Wake 1987). The Bolitoglossa 
subpalmata Species Group (Parra-Olea et al. 2004) is 
an example of this pattern, containing eight species (В. 
bramei, B. gomezi, B. gracilis, B. kamuk, B. pesrubra, B. 
splendida, B. subpalmata, and B. tica) that are distributed 
in the highlands of Costa Rica and western Panama 
(Fig. 1), with high species turnover along latitudinal 
and elevational gradients (AmphibiaWeb 2023; Boza- 
Oviedo et al. 2012; García-París et al. 2008). Although 
the phylogenetic relationships of the Bolitoglossa 
subpalmata S.G. have been relatively well-studied 
(Boza-Oviedo et al. 2012), there 15 still little information 
on the distribution ranges of the species in this group, 
especially B. pesrubra, and the area between the páramo 
of the Cerro Chirripó and the páramo of the Cerro Kamuk 
has not been sampled at all. 

During fieldwork in the subalpine rain paramo in the 
Cordillera de Talamanca, specifically the summits of the 
peaks Durika, Arbolado, Haku, and Utyum (Fig. 1), we 
found some specimens of the B. subpalmata S.G., filling 
in the distribution range for the group. Based on the 
molecular and morphological data reported herein, we 
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describe these samples as a new species of Bolitoglossa 
closely related to B. kamuk and B. pesrubra of the B. 
subpalmata S.G. 


Materials and Methods 


Taxon sampling. In July 2013, March 2015, January 
2016, and October 2016, we carried out collecting trips to 
Cerro Utyum (July 2013 and March 2015), Cerro Dürika 
(January 2016), Cerro Arbolado (October 2016), and 
Cerro Hakü (October 2016) on the Talamanca Mountain 
range within La Amistad International Park, Costa Rica 
(Fig. 1). These trips involved walking a transect from 
Olán, Buenos Aires (9.2788? N, 83.2156? W, 1,500 m 
asl; all GPS coordinates in WGS84 datum) on the Pacific 
slope to the summits of the peaks. The salamanders 
were found by opening bromeliads from previously 
cut-off trees, and a total of 19 specimens referred to 
Bolitoglossa subpalmata S.G. were found in the sampled 
peak summits. 

All the specimens collected for this study were 
humanely euthanized using a topical anesthetic. 
Following the euthanization process, a small tissue 
sample, either liver or the tip of the tail, was taken and 
stored in 96% ethanol. The specimens were fixed in a 
10% formalin solution and transferred to 70% ethanol 
for long-term preservation. Individuals were deposited 
in the herpetological collection of Museo de Zoología 
at Universidad de Costa Rica (UCR) (Appendix 1). 
Museum collection acronyms follow Frost (2023), 
with the additions of EAP (Erick Arias field numbers), 
CRARC (Costa Rica Amphibian Research Center private 
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collection), and NV= No voucher. 


Amplification and sequencing. Total genomic DNA 
was extracted from the ethanol-preserved tissues of 24 
Bolitoglossa specimens using the phenol-chloroform 
standard protocol (Sambrook and Russell 2006). The 
large subunit ribosomal RNA (16S) and cytochrome b 
(cyt Б) mitochondrial genes were amplified. The primers 
16Sar and 16Sbr (Palumbi et al. 1991) were used for 16S, 
and primers MVZ15 and MVZ16 (Moritz et al. 1992) 
were used for cyt b. PCR amplifications were performed 
using a total volume of 15 uL, which contained 1 uL 
DNA template (at 50 ng uL), 0.75 U Taq polymerase 
(Amplificasa®, Biotecnologias Moleculares), 1X PCR 
buffer with 1.5 mM MgCl, 0.2 mM deoxynucleotide 
triphosphates (dNTPs), and 0.3 uM forward and reverse 
primers. The PCR conditions were as follows: for 16S, 
an initial cycle of 5 min at 94 °C, followed by 35 cycles 
of 45 s at 94 °C, 30 s at 55 °C, 45 s at 72 °C, plus a 
final step of 3 min at 72 °С; and for cyt b, an initial 
cycle of 2 min at 94 °C, followed by 38 cycles of 30 s 
at 94 °C, 1 min at 48 °C, 1 min at 72 °C, plus a final 
step of 8 min at 72 °C. The PCR products were cleaned 
with ExoSap-IT (USB Corporation) and sequenced in 
both directions using the original amplification primers 
and BigDye termination reaction chemistry (Applied 
Biosystems). The cycle-sequencing products were 
column-purified with Sephadex G-50 (GE Healthcare) 
and run on an ABI 3500xL Genetic Analyzer (Applied 
Biosystems). Consensus sequences for each individual 
were constructed using SEQUENCHER 5.3 (Genes 
Codes Corp.). The resulting sequences were deposited in 
GenBank (Appendix 1). 


Phylogenetic analyses. The sequences obtained here 
were compared with the sequences of the 16S and cyt b 
mitochondrial genes for 78 specimens of the Bolitoglossa 
subpalmata species group. Sequences of B. aurae were 
used as outgroup and those from B. compacta were used 
to root all trees based on the results in Rovito et al. (2015). 
The list of vouchers and GenBank accession numbers 
used in this study are provided in Appendix 1. Sequence 
alignments were performed using the MUSCLE 3.7 
software (Edgar 2004) with default parameters and 
trimmed to the point where most of the taxa had sequence 
data. PartitionFinder v2.1.1 software (Lanfear et al. 
2017) and the Bayesian Information Criterion (BIC) 
were used to select the best partition scheme and the 
best model of sequence evolution for each partition. A 
single set of branchlengths was used across all partitions 
(branchlengths-linked), and the search for the best 
partition scheme used a heuristic search (scheme=greedy, 
Lanfear et al. 2012). Four subsets were defined a priori: 
one for 16S and three for cyt 5 (partitioned by codon 
position). 

Phylogenetic analyses were performed using both 
Maximum Likelihood (ML) and Bayesian Inference 
(BI) methods. The maximum likelihood analysis was 
performed using Garli 2.01 (Zwickl 2006). To find 
the best tree, ten search replicates were run with the 
following default setting values: streefname — random, 
attachmentsper-taxon = 24, genthreshfortopoterm = 
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100,000, and significanttopochange = 0.00001. For 
bootstrapping, 1,000 pseudoreplicates were run with 
the previous settings and with the following changes: 
genthreshfortopoterm = 10,000, significanttopochange 
=0.01, and treerejectionthreshold = 20, as suggested in 
the Garli manual to speed up the bootstrapping. From 
these bootstraps, a majority rule consensus tree was 
obtained using Sumtrees (Sukumaran and Holder 2010a) 
from DendroPy package version 4.4.0 (Sukumaran and 
Holder 2010b). Bayesian phylogenetic analysis was 
performed using MrBayes 3.2.6 (Ronquist et al. 2012) 
with the partition scheme and the model of sequence 
evolution for each partition as selected previously. Two 
separate analyses were run, each consisting of 50 million 
generations, sampled every 1,000 generations, and four 
chains with default heating parameters. A time-series plot 
of the likelihood scores of the cold chain was examined 
to check stationarity using Tracer 1.6 software (Rambaut 
et al. 2014). The first 25% of trees were discarded as 
burn-in and the remaining trees were used to estimate 
the allcompat consensus tree along with the posterior 
probabilities for each node and each parameter. The Garli 
and MrBayes analyses were run on the CIPRES portal 
(Miller et al. 2010). Genetic distances (uncorrected 
p-distances) were computed using MEGA6 (Tamura et 
al. 2013). 


Morphometrics. The measurements used herein follow 
those proposed by Kubicki et al. (2022). The specimens 
from new localities were compared morphometrically 
only with the species B. kamuk and B. pesrubra, given 
that these are the species of the group that inhabit the 
subalpine paramo and occur geographically close 
(Fig. 1, Appendix 2). The following 29 morphological 
measurements were taken on each of the 31 adult 
specimens of Bolitoglossa subpalmata S.G. that were 
examined in this study: standard length (SL), shoulder 
width (ShW), head width (HeW), neck width (NeW), 
eye width (EW), snout length (SnL), jaw to snout length 
(JSL), lateral gular fold to tip of snout (LGFS), internarial 
distance (IND), naris to lip distance (NLP), intercanthal 
distance (ICD), hind limb length (HLL), front limb length 
(FLL), trunk width (TW), midventral gular fold to snout 
length (VGS), front limb to snout distance (FSL), ulna 
and hand length (UHL), axilla to groin length (AGL), 
vent length (VL), hand width (HaW), hand length (HaL), 
length of Finger III (LF3), width of Finger III (WF3), 
length of Finger II (LF2), foot width (FoW), foot length 
(FoL), length of Toe III (LT3), width of Toe III (WT3), 
and length of Toe II (LT2). The additional proportions 
reported here include: IND/HeW, HeW/AGL, SnL/HeW, 
HaL/VGS, FoL/VGS, HaW/HeW, FoW/HeW, LT2/FoL, 
LF2/HaL, WT3/FoW, WF3/HaW, HaL/HaW, and FoL/ 
FoW. 

The measurements were taken with a Mitutoyo 
Absolute 4084 digital caliper, and rounded to the nearest 
0.1 mm. Given that there are very significant overlaps in 
the numbers of teeth found in the different Bolitoglossa 
species known to occur in Costa Rica (Boza-Oviedo et al. 
2012; García-París et al. 2008; Savage 2002), the number 
of teeth in the specimens examined in this study were not 
counted. We believe that attempting to accurately count 
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the number of teeth in such small specimens (most of 
which were fixed with their mouths in a closed position), 
requires a high-level of invasive manipulation that would 
most likely cause significant and unnecessary damage 
to the mouth cavity. The limb interval is equal to the 
number of costal folds between the tips of the longest 
digits of the adpressed front and hind limbs, expressed in 
0.5 increments (e.g., 4 or 4.5). 


Morphometric statistics. To avoid allometric effects 
relative to the differences in the sizes and shapes 
between species and between individuals, the data 
were transformed using the method of Lleonart et al. 
(2000). In this method, a logarithmic transformation 
of the continuous variables is performed to reduce the 
extreme values. All transformed variables were used in 
the allometric transformation by means of the equation: 


Ж! 


where Y* corresponds to the value of each ofthe dependent 
variables corrected for size and shape, Y, corresponds to 
the value of each dependent morphometric variable, X, is 
the average of the SL variable for all populations, X, is the 
SL value for each individual, and Ё is the regression line 
intercept with the Y-axis resulting from the regression of 
each dependent variable with X,. The intercept is used 
as an allometric transformation factor and is unique for 
each variable. 

A discriminant analysis was performed to determine 
whether the morphometric variables were effective in 
predicting the species. Only those variables that were 
different among at least two of the three putative species 
analyzed were included in the discriminant analysis. The 
following variables were finally used in the discriminant 
analysis: HeW, JSL, ICD, HLL, VGS, HaW, HaL, HaL/ 
VGS, FoL/VGS, and HaW/HeW. The discriminant 
analyses were separated by sex, and were performed 
using R v3.3.3 (R Development Core Team 2013). 


Results 


Molecular analyses. The mitochondrial dataset includes 
samples of all described species in the Bolitoglossa 
subpalmata S.G. The resulting data matrix had a total 
of 80 sequences (31 salamanders collected in this study 
and 49 sequences downloaded from Genbank) with a 
sequence length of 1,334 bp, including gaps, with 527 bp 
for 16S and 807 for cyt b. Three partitions were identified 
with the following substitution models: HK Y+I+G for 
16S + cyt b codon position 3, HKY- for codon position 
1 of cyt b, and TRN+G for cyt b codon position 2. 

The phylogenies from Garli and MrBayes were 
relatively discordant in the internal topology within 
the B. subpalmata S.G. (Fig. 2). However, all 
phylogenetic analyses found the B. subpalmata S.G.to be 
monophyletic, with three internal clades relatively well- 
supported. The first clade was formed by B. bramei and 
B. gomezi, two species that are restricted to premontane 
and montane forests of the Pacific slope. The second 
clade was formed by B. gracilis sister to B. tica + B. 
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subpalmata, and these species are distributed in the 
Volcanic Central Mountain range and the northern part 
of the Cordillera de Talamanca. Finally, a well-supported 
clade was formed by B. kamuk sister to B. pesrubra + 
unnamed taxon, and these three species are restricted to 
the subalpine rain paramo of the Cordillera de Talamanca. 
The phylogenetic position of B. splendida is uncertain. In 
the Bayesian analysis, this species was weakly supported 
as sister to the clade formed by the three species of the 
paramo. In the ML analysis, B. splendida was weakly 
supported as sister to the clade formed by B. bramei + B. 
gomezi, forming the sister clade to the clade formed by B. 
gracilis, B. subpalmata, and B. tica. 

The mitochondrial genetic distances are shown in 
Table 1. Genetic distances between the specimens of the 
unnamed taxon from the subalpine rain paramo and all 
other members of the B. subpalmata species group are 
1.45—4.71% for 16S and 5.03—9.09% for cyt b. 


Morphometric analyses. The specimens examined 
included 13 specimens (8? and 54) of the unnamed 
taxon from the subalpine rain paramo, 15 specimens of B. 
pesrubra (79 and 8), and three specimens of B. kamuk 
(19 and 245). The morphometric variation between the 
three species that inhabit the subalpine rain paramo of 
Costa Rica is shown in Table 2. The discriminant analysis 
correctly classified 100% of the specimens to the species 
(Fig. 3) in both sexes, showing a clear separation between 
the specimens of the new species described below and 
the specimens of B. kamuk and B. pesrubra. The ratios 
HaL/VGS, HaW/HeW, and IND/HeW differed among 
the three species of the subalpine rain páramo. 


Description of New Species 


Bolitoglossa bolanosi sp. nov. 
Bolafios’ Web-footed Salamander 
(Figs. 4—6) 


urn:lsid:zoobank.org:act:3F7E0D7E-906F-4FCD-B4CE-0A48F134DEA1 


Holotype. UCR 22965, an adult male from Costa Rica: 
Provincia de Puntarenas: Cantón de Buenos Aires: 
Distrito de Buenos Aires: the summit of Cerro Arbolado, 
Parque Internacional La Amistad, (9.320°, -83.216°; 
2,600 m asl), collected by Erick Arias and Omar Zúñiga 
on 19 October 2016. 


Paratopotype. UCR 22964, a subadult male, same data 
as holotype. 


Paratypes. UCR 22424, an adult male; UCR 22423, 
an adult female; UCR 22425, a subadult female; UCR 
22422 and UCR 22426, subadult males; and UCR 22427, 
a juvenile from Costa Rica: Provincia de Puntarenas: 
Canton de Buenos Aires: Distrito de Buenos Aires: the 
summit of Cerro Utyum, Parque Internacional LaAmistad, 
(9.323°, -83.187°; 2,870 m asl), collected by Erick Arias, 
Gerardo Chaves, Olmer Cordero, and Omar Zúñiga 
on 30 March 2015. UCR 22421, an adult female from 
Costa Rica: Provincia de Limón: Cantón de Talamanca: 
Distrito de Telire: the summit of Cerro Utyum, Parque 
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Fig. 2. Bayesian phylogenetic inference shows the relationships of the Bolitoglossa subpalmata species group based on the 16S 


and cyt b mitochondrial DNA gene fragments. Bootstrap proportion 
numbers below the branches are posterior probabilities (multiplied 


s from maximum likelihood are shown above the branches, and 
by 100) from the MrBayes analysis. The scale bar refers to the 


estimated substitutions per site. The support values of any node within the species are not shown. The asterisks represent support 


>99. NV = no voucher. 


Internacional La Amistad, (9.333°, -83.180°; 2,913 m 
asl), collected by Erick Arias, Gerardo Chaves, Olmer 
Cordero, and Omar Zúñiga оп 30 March 2015. UCR 
22745, an adult male; UCR 22741-4, adult females; and 
UCR 22746, a juvenile from Costa Rica: Provincia de 
Puntarenas: Cantón de Buenos Aires: Distrito de Buenos 
Aires: the summit of Cerro Haku, Parque Internacional 
La Amistad, (9.322°, -83.203*; 2,660 m asl), collected 
by Erick Arias and Omar Zúñiga on 28 December 2015. 
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UCR 24245, an adult female; UCR 24246, an adult 
male; UCR 24247, a subadult male; and UCR 24248, 
a juvenile from Costa Rica: Provincia de Puntarenas: 
Cantón de Buenos Aires: Distrito de Buenos Aires: 
the summit of Cerro Dürika, Parque Internacional La 
Amistad, (9.3749, -83.303°; 3,240 m asl), collected by 
Omar Zúñiga on 13 January 2016. 


Generic Placement. Assigned to the genus Bolitoglossa 
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= = Bolitoglossa bolanosi differs from its closest relatives 
e = as follows: Bolitoglossa kamuk Boza-Oviedo et al., 2012 
с А Š = is smaller and more slender with mean SL 34.6—38.4 mm 
5 CE сы (larger and robust, SL = 39.43-50.01 тт); internarial 
= cs = eb & a M PONE distance relatively shorter, IND/HeW = 0.26+0.05 
Ф ОС TE E EC ee ee (internarial distance relatively longer, IND/HeW = 
Е: cj cá ч c с c cd c с 0.31+0.05); shorter tail TL/SL 98-99% (tail long, TL/SL 
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Table 2. Morphometric data for Bolitoglossa bolanosi sp. nov., B. kamuk, and B. pesrubra. These data were obtained prior to 
morphometric correction. 


B. bolanosi sp. nov. B. kamuk 
LEBER MeanzS.D. (Range) Meanz S.D. (Range) MeansS.D. (Range) 
[... | 9e $(N-5 4-2 | 9(WN-) S (N=8) 


43.05+5.54 41.83+6.55 35.2+0.85 50.6+3 48.34+3.61 
(32.11-50.01) | (34.01-49.04) (34.6-35.8) (47.6—54.3) (47.6—54.78) 
SE 5.64+0.98 5.34+1.02 6.3+0.69 5.46+0.52 
(4.25-7.49) (4.4-7.09) (3.8-4) (5.2-7.3) (5.2-6.43) 
mes 6.74+0.7 6.43+0.72 5.1540.21 7.17+0.75 7.02+0.41 
(5.45-7.7) (5.46-7.3) (5-5.3) (6.09-8.3) (6.09-7.76) 
Tos 5.62+0.84 5.17+0.74 4.05+0.07 6.2+0.77 5.38+0.52 
(3.99-6.72) (4.26-6.17) (44.1) (5.3-7.4) (5.3-6.53) 
= 2.32+0.28 2.5+0.48 2.15+0.21 2.56+0.22 2.5140.22 
(1.78-2.59) (2.01-3.18) (0953) (22:908) (222218) 
a 2.75+0.32 2.88+0.43 
(2.26-3.09) (2315343) 
m 740.64 740.88 
(5.99-8.15) (6.01-8.21) 
10.99+1.19 10.75+1.56 
LOES (8.46-12.17) (9.06—12.9) 
1.94+0.14 2.25+0.56 
(1.74-2.12) (1157511) 


B. pesrubra 


3.9+0.14 


— 
N 


10.30.28 
(10.1-10.5) 


= ] 
8.4 
3 1.95+0.07 3.03+0.18 30.76 
(1.9-2) (2.8-3.3) (2.8-3.78) 
5.95+0.49 8.09+0.5 8.08+0.48 
(5.6-6.3) (7.6-9) (7.6—8.83) 
8.95+0.21 11.84+0.65 11.94+0.63 
18 2:2 72:0: 19 2.61+0.36 
| (2—2.5) (2—3.16) 
NLP 0.9+0.09 0.87+0.15 1 0.75+0.07 1.07+0.13 1.05+0.16 
(0.74—1.01) (0.74—1.11) (0.7—0.8) (0.9-1.3) (0.9—1.28) 
ICD 3.3+0.34 3.27+0.38 28 2.45+0.07 3.14+0.57 3.51<0.28 
(2.88—3.8) (2.93-3.88) j (2.4-2.5) (1.9-3.6) (1.9-3.98) 
HLL 111.62 11.07+1.54 04 8.3+0.28 12.26+0.85 12.711.36 
(8.03-12.79) (9.47—13.09) | (8.1—8.5) (11.2—13.6) (11.2—14.16) 
TW 6.62+1.39 5.26+0.63 54 3.9+0.14 6.93+0.96 6.1+0.46 
(5.16—9.2) (4.3—6.02) j (3.8—4) (5.8—8.4) (5.8-6.9) 
10.38+1.54 11.28+3.87 
(7.75-11.91) (8.06—17.62) E N 217 (10.1-12.4) (10.1—13.61) 
(11-122) 
FSL 12.75+1.42 12:59:51:77 19 14.46+0.79 14.37+0.85 
(10.16-14.8) (10.91—15.17) i (13.5—15.6) (13.5—15.59) 
UHL 7.2242 6.5122:65 6.5 8.1+0.64 6.16+2 
(2.63-8.58) (3.29.01) (7.4—9.4) (7.4-9.2) 
AGL 21.94+3.75 2102224109 93 17.9+0.28 27 371580 24.71+1.96 
(16.14—26.4) (16.18—24.13) (17.7-18.1) (25.2-30.1) (25.2-28.35) 
VL 4.34+0.7 4.06+0.85 | 2.56+0.62 4.86+0.53 4.36+0.64 
(3.05-5.09) (2.93—4.94) (2.12-3) (4.4-5.8) (4.4—5.78) 
HaW 3.66+0.5 3.79+0.84 3 2.65<0.21 4.24+0.44 4.36+0.57 
(3.04-4.41) (2.68—4.84) i 
3.85+0.64 3.811 
(2.75-4.64) | (291-535) 
0.86+0.09 0.79+0.15 
(0.73-1.05) | (059-09%) EN 
LF? 1.6140.22 1.8140.29 16 1.45+0.07 2.01+0.22 2.24+0.4 
(1.33-1.92) (1.59-2.31) | (1.4—1.5) (1.8—2.4) (1.8-2.7) 
LF3 2.02<0.25 2.29+0.45 23 П 5021 2.56+0.24 2.76+0.42 
(1.75-2.39) (1.86—3) | (1.6—1.9) (2.4-3.1) (2.4-3.45) 
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(8.8-9.1) (112-129) (11.2-12.96) 
11.2740.74 [1515-1-22 
ОС 10.49+1.13 10.47+1.81 8.45+0.35 11.66+0.52 11.40.56 
(8.18-11.68) (872132) (82-87) (01-10 4) 
(2528) (3855 1) (318251) 
2.6+0.14 4.04+0.43 4.2+0.49 
(25209) (3:855) (3.84.9) 


0.77+0.08 0.82+0.11 
(0.7-0.9) (0.7—0.99) 


(N 
3 
2 
o 
1 
1 
1 
1 
3 
3 
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Table 2 (continued). Morphometric data for Bolitoglossa bolanosi sp. nov., B. kamuk, and B. pesrubra. These data were obtained 
prior to morphometric correction. 


[КЕ | Mean+S.D. (Range) Mean+S.D. (Range) Mean+S.D. (Range) 
| | 208 dN-5 | N=) | &@=2) 9(N=7) d N=8) 


m 4.44+0.76 4.56+0.87 - 3.4540.21 4.96+0.4 5.17+0.77 
(3.44—5.42) (3.63-5.74) | (3.3-3.6) (4.5-5.6) (4.5-6.02) 
4.34+0.58 4.23+0.91 4.6+0.35 4.78+0.39 
(3.42-5.03) (3.27-5.53) (4.1-5.2) (4.1-5.3) 
0.9:0.16 0.77+0.1 0.65+0.07 0.8+0.1 0.86+0.13 
(0.68—1.24) (0.67-0.92) (0.6-0.7) (0.7-1) (ОИ) 
1.8240.19 1.93+0.16 T 1,240.14 2.3140.22 2.4440.5 
(1.52-2.05) (25) | (ТЭ) (222 2) (2-3.09) 
2.44+0.44 2.46+0.22 E: 1.65+0.07 2.79+0.17 3.05+0.42 
(1.66-2.94) (2.15-2.68) | (1.6-1.7) (2.6-3) (2.6-3.73) 
0.24+0.01 0.25+0.01 0.24+0 0.23+0.01 0.24+0.01 
NEL (0.23-0.26) (0.23-0.27) (0.24-0.24) (0.22-0.24) (0.22-0.25) 
0.29+0.02 0.35+0.05 0.32+0.03 0.37+0.04 
а (0.26-0.32) (0.29-0.43) (0.29-0.38) (0.29-0.42) 
0.51+0.03 0.51+0.05 0.54+0.01 0.51+0.02 
QUE (0.47-0.58) (0.46-0.58) pene (0.53-0.55) (0.53-0.54) 
0.16+0.01 0.15+0.01 0.15+0.01 0.14+0.01 0.15+0.01 
HET. (0.15-0.17) (0.15-0.17) 0186 (0.14-0.15) (0.12-0.16) (0.12-0.16) 
0.31+0.03 0.31+0.02 кк 0.29+0.01 0.26+0.02 0.28+0.01 
(0.25—0.36) (0.28—0.34) j (0.28-0.29) (0.23-0.29) (0.23-0.3) 
0.45+0.03 ы 0.38+0.03 0.42+0.03 0.43+0.1 
36-0. (0.41-0.48) | (0.36-0.4) (0.37-0.46) (0.37-0.49) 
0.26+0.01 0.27+0.01 т 0.24+0.01 0.24+0.01 0.26+0.02 
(0.25-0.27) (0.26-0.28) | (0.23-0.24) (0.22-0.26) (0.22-0.3) 
0.24+0.01 0.26+0.05 e 0.23+0.01 0.22+0.02 0.24+0.02 
(0.22-0.26) (0.23-0.26) | (0.22-0.23) (0.19-0.24) (0.19-0.27) 
0.3740.03 0.36+0.03 re 0.3140.03 0.35+0.03 0.3740.03 
(0.33-0.4) (0.33-0.41) (0.29-0.33) (0.31-0.4) (0.31—0.43) 
0.4140.02 0.4+0.03 "m 0.36+0.01 0.39+0.02 0.42+0.03 
(0.37—0.44) (0.37-0.43) | (0.34—0.37) (0.37-0.42) (0.37-0.45) 
0.54+0.04 0.58+0.07 m 0.52+0.06 0.59+0.04 0.62+0.07 
(0.48—0.59) (0.49—0.66) | (0.47-0.56) (0.54—0.66) (0.54—0.76) 
0.66+0.06 0.7+0.06 те 0.67+0.07 0.7+0.07 0.73+0.09 
(0.55-0.73) (0.64-0.79) | (0.62-0.72) (0.63-0.84) (0.63-0.9) 
0.27+0.03 0.31+0.03 035 0.23+0.02 0.33+0.04 0.35+0.07 
(0.24—0.3) (0.26-0.32) (0.222025) (0.27-0.39) (0.27-0.45) 
0.55+0.05 0.56+0.08 m 0.48+0.05 0.56+0.05 0.60.08 
(0.48-0.63) (0.47-0.64) | (0.44—0.52) (0.51-0.65) (0.51-0.71) 
0.42+0.03 0.47+0.09 0.4+0.05 0.5+0.05 0.5101 
DL (0.37-0.48) (0.35-0.54) nt (0.37-0.43) (0.42-0.57) (0.42-0.68) 
0.43+0.08 0.49+0.08 0.56+0.06 0.540.05 0.53+0.08 
LEAD (0.29-0.58) (0.42-0.6) Oe (0.52—0.6) (0.45-0.58) (0.45-0.63) 
0.2+0.03 0.17+0.01 0.19+0.03 0.16+0.03 0.17+0.02 
ЕИ (0.16-0.24) (0.16-0.18) T (0.17-0.21) (0.14-0.22) (0.14-0.21) 
0.24+0.03 0.2140.03 0.23+0.02 0.18+0.01 0.19+0.03 
Ши (0.21-0.28) (0.16-0.26) 2273 (0.21-0.24) (0.17-0.21) (0.17-0.24) 
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102.3-117.2%); dorsal ground color relatively uniform, 
orange to black (dorsal ground usually mottled or blotchy); 
males with rounded snout (snout strongly truncated in 
males). Bolitoglossa pesrubra (Taylor, 1952) has hands 
and feet more webbed (Fig. 5C—D), usually less than the 
ultimate phalange free (with at least the ultimate phalange 
free beyond the interdigital tissue margin, Fig. 5A-B); has 
red color on forelimbs and usually on hind limbs (dorsum 
variable but never with red on hind limbs or forelimbs). 
Hands relatively wider, HaW/HeW = 0.61+0.07 (hands 
narrower, HaW/HeW = 0.56+0.06). 


Description of holotype. Adult male having a SL of 
49.05 mm (Fig. 4). Head slightly wider than neck and 
shoulders (HeW 7.3 mm, NeW 6.2 mm, ShW 7.1 mm), 
with the greatest width of the head just posterior to the 
articulation of the jaws; truncate in dorsal outline and 
rounded to truncate in profile; snout moderate (SnL 3.4 
mm, 6.9% of SL), with nearly terminal non-protruding 
small nostrils (LNH 0.6 mm, RNW 0.4 mm) directed 
anterolaterally; internarial area convex in dorsal outline. 
Snout protruding beyond the anterior margin of the lower 
lip in lateral view. Eyes relatively large (EW = 92% of 
SnL), weakly protruding beyond the dorsal and ventral 
outline of the head, directed anterolaterally, with a 
distinct suborbital groove. Top of head flat and smooth, 
tapering slightly toward anterior terminus, lacking 
contrasting interorbital or other dermal structures. 
Canthus rostralis weakly rounded; intercanthal area flat 
to slightly convex; and loreal region slightly concave. 
Evident cirri (nasolabial protuberances) on tip of snout; 
nasolabial grooves start at ventrolateral margins of nares 
and terminate in a cirrus. Gular fold is well-defined, 
starting on the dorsolateral portion of the neck, below the 
postorbital groove. Evident mental gland is visible under 
the skin of the anterior intermandibular region. 

Arms relatively long and slender (FLL = 12.6 mm, 
25.7% of SL), without noticeable hypertrophied forearm 
compared to the upper arm. Hands well-developed and 
slender (HaL = 5.4 mm, 40.9% of VGS; HaW = 4.8 mm, 
65.8% of HeW). Fingers II, III, and IV protrude freely, 
with at least the ultimate phalange free beyond interdigital 
tissue margin (LF2 2.3 mm, LF3 3.0 mm; Fig. 5). Tips of 
fingers rounded; terminal pads weakly discernible on the 
ventral surfaces of fingers. Relative lengths of fingers on 
right hand: I < IV < II < III. 

Legs moderately long and slender (HLL 13.1 mm, 
26.7% of SL). Feet well-developed and slender (FoL 
5.5 mm, 41.7% of VGS; FoW 5.7 mm, 78.1% of 
HeW). Toes II, III, IV, and V protruding freely beyond 
interdigital tissue margin (LT2 1.9 mm, LT3 2.7 mm), 
toe I with minimal indentation at interdigital spaces. 
Toe III is most free of interdigital tissue, with about the 
entire distal phalanx protruding. Tips of toes rounded; 
terminal pads weakly discernible on ventral distal 
surface of toes. Relative lengths of toes on right foot: I 
МОШ ЕИ 

Body subcylindrical (slightly wider than high) in 
cross-section, and relatively slender (TW = 6.0 mm; TW 
= 26.5% of AGL). Between the axilla and groin, 11 costal 
grooves are visible, 13 if counting axillary and inguinal 
grooves; costal grooves are most visible on ventral and 
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lateral portions of the body. Adpressed limbs separated 
by one costal fold; 12 costal folds total between axilla 
and groin. The tail is long, cylindrical in cross-section, 
with an evident constriction at the base, and some caudal 
grooves discernible on the anterior portion of the tail. The 
skin on the surfaces of the head, body, limbs, and tail is 
smooth. 


Coloration in life. The ground color of dorsal surfaces of 
the head, trunk, tail, hind limbs, and forelimbs is brownish 
black to brownish violet with numerous fine lighter 
patches of chromatophores scattered throughout the 
dorsal surface, especially concentrated on the head. The 
head and hindlimbs are lighter than the trunk. A pair of 
irregular dorsolateral stripes, which run from the tail base 
to the intercanthal area crossing the superior eyelids, are 
formed by yellowish-bronze blotches. The iris is bright 
dark bronze with a dark brownish-black reticulation. 
The upper surfaces of the arms are lighter than the trunk, 
pinkish brown, with orange blotches in the proximal 
portion of the humerus and on the ventrolateral surface 
of the ulna. The upper surface of the legs is similar in 
color to the trunk, uniform brownish black, except in the 
foot which is paler. The dorsal and dorsolateral surfaces 
of the tail are nearly uniform brownish black. 

The ventrolateral surfaces of the body, tail, hindlimbs, 
and forelimbs are lighter than the dorsal surface. The 
gular surfaces are paler than the venter, consisting of 
lighter brown with a paler blotch on the anterior part, 
the mental gland. The ventral surface of the trunk and 
the tail are slightly lighter than the dorsal surface. The 
ventral surfaces of the arms and legs are lighter than 
the dorsal surfaces with irregular orange blotches. The 
palmar and plantar surfaces are light brown with black 
chromatophores. 


Coloration in ethanol. After more than seven years in 
ethanol (70%), the overall coloration of the holotype 
has darkened throughout and contains a principal dark 
brown-blackish tone. 


Measurements (in mm), limb interval, and percentages 
of the holotype. SL 49.05; TL 57.5; ShW 7.1; HeW 7.3; 
NeW 6.2; EW 3.2; SnL 3.4; JSL 8.2, LGFS 12.9 LNH 
0.6; RNW 0.4; IND 3.1; NLP 1.1; ICD 3.9; HLL 13.1; 
FLL 12.6, TW 6.0; VGS 13.2; FSL 15.2; UHL 9.0; AGL 
22.8; VL 4.9; HaW 4.8; HaL 5.4; LF2 2.3; LF3 3.0; WF3 
1.0; FOW 5.7; FoL 5.5; LT2 1.9; LT3 2.7, WT3 0.9. Limb 
interval 1. Measurements in relative percentages: VGS/ 
SL 26.9%; IND/HeW 42.596; AGL/SL 46.5%; HeW/SL 
14.9%; HeW/AGL 32.096; SnL/ HeW 46.6%; LNH/HeW 
8.2%; LNH/SL 1.2%; RNW/HeW 5.5%; RNW/SL 0.8%; 
HLL/SL 26.7%; FLL/SL 25.7%; HaL/VGS 40.9%; FoL/ 
VGS 41.7%; HaW/HeW 65.7%; FoW/HeW 78.1%; LT2/ 
FoL 34,596; LF2/HaL 42.6%; WT3/FoW 15.8%; WF3/ 
HaW 20.8%. 


Noteworthy variation. The female paratypes (UCR 
22421 and UCR 24245) have more rounded snouts, 
and much less defined and protruding cirri or nasolabial 
protuberances, which are sexually dimorphic features. In 
coloration, this species is extremely polymorphic, some 
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Fig. 3. Linear discriminant analysis shows the morphological separation between the new species and the two species close to it. 


individuals are almost uniform black, which is common 
in juveniles (UCR 22247, Fig. 6B). The paratopotype 
(UCR 22964, Fig. 6A) has the pair of dorsolateral stripes 
more vivid and continuous, the blotches in the upper 
eyelid are iridescent green; in addition, it has a band 
of pale orange from posterior head to tail, where it is 
suffused with the red of the tail. The presence of red on 
the tail is common; generally the red color in the proximal 
portion is only on the dorsal surface, but the tail becomes 
completely red in the distal portion. Several specimens 
showed a similar coloration morpho “mottled” (Fig. 6C), 
with ground color brownish black to reddish brown with 
numerous and irregular yellow to red blotches or spots; 
these specimens lack the dorsolateral stripes. The female 
UCR 24245 (Fig. 6D) has a striking coloration pattern 
formed by large yellow blotches about a brownish-black 
ground color; although less evident this specimen has the 
pair of dorsolateral stripes but the blotches are larger and 
suffused in the back. 


Measurements (in mm), limb intervals, and 
percentages of the paratypes. SL 39.4—50.0; ShW 
4.9-7.5; HeW 6.3-7.7, NeW 5.0—6.7; EW 1.8-3.2; SnL 
2.4—3.4; JSL 6.0-8.2; LGFS 10.2-12.9; IND 1.7-3.1; 
NLP 0.7-1.1; ICD 2.9-3.9; HLL 9.9-13.1; FLL 8.9— 
17.6; TW 5.2-9.2; VGS 9.6-13.2; FSL 11.4-15.2; UHL 
6.8—9.0; AGL 18.6—26.4; VL 3.5—5.1; HaW 3.0-4.8; HaL 
3.4—5.4; LF2 1.3-2.3; ГЕЗ 1.83.0; МЕЗ 0.7-1.1; FoW 
3.4—5.7; FoL 3.8-5.5; LT2 1.5-2.1; LT3 1.72.9; WT3 
0.7-1.2. Limb intervals 3—5. Measurements in relative 
percentages: VGS/SL 23.4—26.9%; IND/HeW 25.8- 
42.6%; AGL/SL 46.4—57.6%; HeW/SL 14.5—16.9%; 
HeW/AGL 25.2—35 7%; SnL/HeW 35.7-47.6%; HLL/ 
SL 24.6—26.8%; FLL/SL 22.3—26.1%; HaL/VGS 32.7— 
40.5%; FoL/VGS 37.1—44.4%; HaW/HeW 48.0—66.3%; 
FoW/HeW 54.5—78.6%; LT2/HeW 24.1—31.9%; LT3/ 
FoW 46.7-63.0%; LT2/FoL 34.9-45.8%; LF2/HaL 29.3— 
46.9%; WT3/FoW 15.5—24.1%; WF3/HaW 16.2—28.3%. 


Etymology. The name *bolanosi" is a patronym honoring 
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the Costa Rican herpetologist Federico Bolafios, and is 
used as a noun in the genitive case. We name this species 
after our dear friend in recognition of his scientific 
contributions to the knowledge of the herpetology of 
Costa Rica, as curator of the Herpetology section at 
Museo de Zoología of Universidad de Costa Rica, and as 
the mentor of most herpetologists present in the country. 


Habitat and natural history observations. The habitat 
of Bolitoglossa bolanosi in the subalpine rain páramo is 
characterized by having a very short dry season (one to 
two months), annual precipitation ranging from 1,000 
to 2,000 mm, and annual temperatures between 3 and 6 
°C. The paramo vegetation consists of extensive, deep 
moss mats, spongy soil, ferns, and small isolated trees 
with arboreal bromeliads. Bolitoglossa bolanosi was 
found within both moss and bromeliads. It also occurs 
in Montane rain forest which is dominated by Quercus 
covered with moss, with a great abundance of bryophytes 
and epiphytes. The type locality occurs in the Montane 
rain forest (Bolafios et al. 2005; Holdridge 1967), 
characterized by a very short dry season (one to two 
months), an annual precipitation range of 2,200 to 4,500 
mm, and annual temperatures from 6 to 12 °C. 

Very little is known about the natural history of B. 
bolanosi, but it is important to note that females guarding 
a clutch of eggs were found on Cerro Hakü (December 
2015) and Cerro Durika (January 2016), both under moss 
on the floor. All specimens from Cerro Utyum were found 
within bromeliads at heights of >3 m over the soil, despite 
active searches in the moss of the páramo. Unlike in the 
peaks Dürika, Arbolado, and Haku, the specimens were 
mainly found in the moss at «2 m over soil, despite active 
searches in the bromeliads at >3 m above the ground 
(especially on Cerro Arbolado and Cerro Haku). In the 
summit of Cerro Dürika, B. bolanosi is very near (~4 km) 
to a site with B. pesrubra. In the summit of Cerro Utyum, 
it is to only ~7 km to a site with B. kamuk, and in this peak 
B. bolanosi is sympatric with an unnamed miniaturized 
salamander related to B. pygmaea. In addition, the type 
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Fig. 4. Holotype of Bolitoglossa bolanosi sp. nov. (UCR 22965) in life on a white background. Photograph by Erick Arias. 


| 


> E E 


Fig. 5. Dorsal views of the hand and foot of the holotypes of (А-В) Bolitoglossa bolanosi sp. nov. (UCR 22965) апа (C-D) B. 
pesrubra in preservation. Photograph by Erick Arias. 
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Fig. 6. In-life photographs of variations within Bolitoglossa bolanosi sp. nov. (A) Paratopotype male subadult UCR 22966, (B) 
paratype juvenile UCR 24247, (C) paratype adult female UCR 22421, and (D) paratype adult female UCR 24245. Photographs A 


and C by Erick Arias, B and D by Omar Zúñiga. 


locality of B. splendida is only ~4 km from the site with 
B. bolanosi. 


Distribution. The known distribution area of Bolitoglossa 
bolanosi is very small, restricted throughout to ~15 km 
on the summits of the peaks Dürika, Arbolado, Hakü, 
and Utyum on the Cordillera de Talamanca (Fig. 1). The 
altitudinal range of the new species is 2,550-3,240 m 
asl. All the populations of the new species were found in 
primary vegetation (páramo and forest) and all are within 
a protected area, La Amistad International Park. 


Conservation status. The conservation status of this 
species is uncertain; however, its known distribution 
range is small («20,000 km?) and all know populations 
are restricted to summits of peaks, sites that are highly 
vulnerable to forest fires and other effects of climate 
change. We suggest that it should be tentatively 
considered as Vulnerable (VU) following the IUCN 
(2017) category criteria. 


Discussion 


Costa Rica hosts 56 species of salamanders, and it is the 
sixth most diverse country in the World in terms of total 
species richness and the most diverse relative to its total 
area (AmphibiaWeb 2023; Boza-Oviedo et al. 2012; Frost 
2023). Several new species of salamanders have been 
described from ICA in this century (Arias and Kubicki 
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2018; Bolafios and Wake 2009; Boza-Oviedo et al. 2012; 
García-París et al. 2008; Hanken et al. 2005; Kubicki 
2016; Kubicki and Arias 2016; Kubicki et al. 2022; 
Wake et al. 2007). Nevertheless, the taxonomic studies 
of salamanders in this region have been obscured due 
to a conservative morphology (Arias and Kubicki 2018; 
Boza-Oviedo et al. 2012; Kubicki et al. 2022). Recently, 
Kubicki et al. (2022) found low genetic distances among 
species of Nototriton from Costa Rica; however, they 
showed that two different species of Nototriton coexist 
(sympatry) in two sites, providing strong evidence of 
speciation. We also found relatively low (Table 1) genetic 
distances among the species of the B. subpalmata S.G., 
however, we found two sites separated by less than 7 km 
in a straight line from where two different species occur. In 
the first case, in the páramo of the Cerro Dürika, samples 
of B. bolanosi (UCR 24245-7) were found only 4 km 
in a straight line from specimens of B. pesrubra (UCR 
23049). The area separating these two sites consisted of 
homogeneous páramo, lacking any significant features 
that could be viewed as a barrier that would prevent the 
potential mixture of individuals or genes. However, the 
genetic distance among these organisms is 2.14% for 
16S, providing evidence of genetic isolation. A second 
case was found in the páramo of the Cerro Utyum, where 
the B. bolanosi specimens (UCR 22421-7) were found 
only 7 km in a straight line from B. kamuk (UCR 23047). 
Similarly, these sites are separated by homogeneous 
paramo and we found genetic distances of 1.35—1.8% for 
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16S. We acknowledge the limitations of our conclusions, 
since our genetic analyses were based solely on the 
mitochondrial results. However, we think that all the 
evidence supports our hypothesis about the relationships 
within the B. subpalmata S.G. 

Bolitoglossa bolanosi is the fourth species of 
salamander endemic to the páramo of ICA and the 
Montane rain forest surrounding it. The other three 
species also belong to the B. subpalmata S.G. Although 
there are relatively few species, the role of the páramos 
of ICA in the speciation of salamanders is interesting, 
especially since the paramos correspond only to 0.3% 
of the continental surface of Costa Rica (Kappelle and 
Horn 2016). The current isolated patches of páramos in 
ICA and its contraction-expansion dynamics in the past 
can explain (in part) the species formation in this habitat. 
In addition, it is noteworthy that the larger páramo in 
Costa Rica, the Cerro Chirripó, has not been sampled for 
salamanders. This páramo corresponds to two-thirds of all 
páramo in the country. There is a need for more studies in 
the páramos of ICA, especially due to the vulnerability of 
this ecosystem to climate change (Enquist 2002; Nogués- 
Bravo et al. 2007). In the past, B. pesrubra was relatively 
common (Bolafios and Wake 2009; Vial 1968). However, 
its populations have been greatly reduced (Boza-Oviedo 
et al. 2012; García-Paris et al. 2008), and we do not fully 
understand the cause of this reduction. 
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Appendix 1. Institutional voucher numbers, locality information, and GenBank accession numbers for the specimens 
used in the molecular phylogenetic analyses. Museum collection acronyms follow Frost (2023) with the addition of 
EAP to refer to Erick Arias field numbers and CRARC to refer to the Costa Rica Amphibian Research Center private 
collection. NV= No voucher. CR = Costa Rica, HN = Honduras, MX = Mexico. 


B. mexicana MVZI176838 Chiapas, MX 320 17.382 :91:793 GU725457 | GU725470 


? 
T 


B. aurae UCR22842 Turrialba, Cartago, CR 1300 KX779527 | KX779528 


B. gracilis MVZ229170 Paraíso, Cartago, CR 1680 9.692 -83.782 AY526121 AF212067 


B. gracilis MVZ229171 Paraíso, Cartago, CR 1680 9.692 -83.782 AY526122 AF212068 
B. gracilis UCR20541 Paraíso, Cartago, CR 1400 9.734 -83.764 OR162558 


B. tica UCR12065 Guarco, Cartago, CR 2440 9.686 -83.894 AY 526137 AF212089 
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Appendix 1 (continued). Institutional voucher numbers, locality information, and GenBank accession numbers for 
the specimens used in the molecular phylogenetic analyses. Museum collection acronyms follow Frost (2023) with 
the addition of EAP to refer to Erick Arias field numbers and CRARC to refer to the Costa Rica Amphibian Research 
Center private collection. NV= No voucher. CR = Costa Rica, HN = Honduras, MX = Mexico. 


3 — - T T - — ' "| 


ed — —-- 


2300 957.12 -83.858 JQ899162 JQ899192 
1850 10.133 AF212093 


1500 10.033 -84.483 AF212090 
1500 10.033 -84.483 EU448107 AF212091 


1700 9.998 -83.964 NON AF212095 
1500 AF416697 AF212094 


ND ND AF212092 


ND 
1600 10.177 -84.391 ОК162562 CE 
EE OR162561 OR147301 


B. kamuk UCR20852 Talamanca, Limon, CR 3126 D -83.059 10899143 10899173 
B. kamuk UCR20853 Talamanca, Limón, CR 9:255 -83.059 JQ899144 JQ899174 


B. kamuk CR20854 Talamanca, Limón, CR 3126 2:255 -83.059 JQ899145 JQ899175 
Talamanca, Limón, CR 2870 9277 -83.135 OR162582 pu 


S. 
ME. 


| 


B. tica UCR20514 Paraíso, Cartago, CR 


B. subpalmata | MVZ190875 Barva, Heredia, CR 


B. subpalmata | MVZ194827 San Ramón, Alajuela, CR 


B. subpalmata | MVZ194828 San Ramón, Alajuela, CR 


Vazquez de Coronado, San 


MVZ194889 José, CR 


B. subpalmata 


B. subpalmata | MVZ229172 Monte Verde, Puntarenas, CR 


| 
| 
І 


NV-Braulio 


Carrillo Heredia, CR 


B. subpalmata 


B. subpalmata | UCR20833 Alfaro Ruiz, Alajuela, CR 


B. subpalmata | UCR21783 


= 


B. kamuk 


ре 


СК23047 


B. bolanosi 
Sp. nov. 


CR22421 


B. bolanosi 
Sp. nov. 


UCR22422 Buenos Aires, Puntarenas, CR 


B. bolanosi 
sp. nov. 


UCR22423 Buenos Aires, Puntarenas, CR 


B. bolanosi 
sp. nov. 


UCR22424 Buenos Aires, Puntarenas, CR 


B. bolanosi 
Sp. nov. 


CR22425 Buenos Aires, Puntarenas, CR 


B. bolanosi 


UCR22426 Buenos Aires, Puntarenas, CR 
sp. nov. 


B. bolanosi 
sp. nov. 


JCR22427 Buenos Aires, Puntarenas, CR 


B. bolanosi 
Sp. nov. 


JCR22741 Buenos Aires, Puntarenas, CR 


B. bolanosi 
sp. nov. 


UCR22742 Buenos Aires, Puntarenas, CR 


B. bolanosi 
sp. nov. 


UCR22743 Buenos Aires, Puntarenas, CR 


B. bolanosi 
sp. nov. 


UCR22744 Buenos Aires, Puntarenas, CR 


B. bolanosi 
sp. nov. 


UCR22745 Buenos Aires, Puntarenas, CR 


B. bolanosi 
Sp. nov. 


UCR22746 Buenos Aires, Puntarenas, CR 


2660 9:322 -83.203 OR162581 
2600 -83.216 OR162566 
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Appendix 1 (continued). Institutional voucher numbers, locality information, and GenBank accession numbers for 
the specimens used in the molecular phylogenetic analyses. Museum collection acronyms follow Frost (2023) with 
the addition of EAP to refer to Erick Arias field numbers and CRARC to refer to the Costa Rica Amphibian Research 
Center private collection. NV= No voucher. CR = Costa Rica, HN = Honduras, MX = Mexico. 


И — J 


UCR24245 Buenos Aires, Puntarenas, CR 3240 9.374 -83.303 OR162568 OR147304 


[ 


B. bolanosi 
sp. nov. 


B. bolanosi 
sp. nov. 


UCR24246 


B. pesrubra 


CRARCO190 Dota, San José, CR 2750 9.652 -83.849 OR162585 OR147314 
CRARCO192 Dota, San José, CR 2750 9.652 -83.849 OR162586 OR147315 
ND ND 


B. pesrubra 


B. pesrubra 


ND ND 
DBW5118 ND D 


EAP0771 Buenos Aires, Puntarenas, CR 


B. pesrubra 


B. pesrubra MVZ190928 2760 | 906 | s2 | | AF212072 
СИСТИНА ЕСЕСІНЕ 
eae [мами mm [om m pem [ - pe 
элин [маны [emet | ж [эш [9m | mm 


{ш н лу уз: 2980 9.617 -83.833 EE AF212077 


Me mm ШИШЕНИН T E [eee [mme 


B. pesrubra 


B. pesrubra 


Buenos Aires, Puntarenas, CR 3240 9.374 -83.303 OR162569 OR147305 


B. pesrubra 
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B. pesrubra 
B. pesrubra 
B. pesrubra 
B. pesrubra 


B. pesrubra 
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Appendix 2. Specimens examined in the morphological analysis. All voucher numbers below are ‘UCR’ 
numbers, and the specimens were housed at the Museo de Zoologia at the Universidad de Costa Rica. 


Bolitoglossa bolanosi sp. nov. 


COSTA RICA: Limón, Talamanca, Cerro Utyum: (22421 9); Puntarenas, Buenos Aires, Cerro 
Arbolado: (229655); Puntarenas, Buenos Aires, Cerro Наки: (242459, 242464); Puntarenas, 
Buenos Aires, Cerro Haku: (22741—49, 227450); Puntarenas, Buenos Aires, Cerro Utyum: 
(2242245, 224239, 224243, 224259). 


Bolitoglossa kamuk 


COSTA RICA: Limón, Talamanca, Cerro Diká: (230479); Limón, Talamanca, Cerro Kamuk: 
(20852-38). 


Bolitoglossa pesrubra 


COSTA RICA: Cartago, El Guarco, Cerro de la Muerte: (239124, 224369); Cartago, Paraíso, Cerro 
de la Muerte: (138104, 190344, 191259, 190999, 236344); San José, Dota, Cerro de la Muerte: 
(120674, 207939, 20797-99); San José. Pérez Zeledón, Cuericí: (117124, 117574, 121098). 
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